Spin squeezing and pairwise entanglement for symmetric multiqubit states 
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We show that spin squeezing implies pairwise entanglement for arbitrary symmetric multiqubit 
states. If the squeezing parameter is less than or equal to 1, we demonstrate a quantitative relation 
between the squeezing parameter and the concurrence for the even and odd states. We prove that 
the even states generated from the initial state with all qubits being spin down, via the one-axis 
twisting Hamiltonian, are spin squeezed if and only if they are pairwise entangled. For the states 
generated via the one-axis twisting Hamiltonian with an external transverse field for any number 
of qubits greater than 1 or via the two-axis counter-twisting Hamiltonian for any even number 
of qubits, the numerical results suggest that such states are spin squeezed if and only if they are 
pairwise entangled. 
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INTRODUCTION 



MM 

correlated 

states with reduced fluctuations in one of the collective 
spin components, with possible applications in atomic 
interferometers and high precision atomic clocks. It is 
found that spin squeezing is closely related to and implies 
quantum entanglement 21, 22, ,i23| . As there are various 
kinds of entanglement, a question naturally arises: what 
kind of entanglement does spin squeezing imply? Re- 
cently, it has been found that, for a two-qubit symmetric 
state, spin squeezing is equivalent to its bipartite entan- 
glement |23|; i.e., spin squeezing implies bipartite entan- 
glement and vice versa. Here, we generalize the above 
result to the multiqubit case, and study relationships be- 
tween spin squeezing and quantum entanglement. 

Specifically, we first show that spin squeezing implies 
pairwise entanglement for arbitrary symmetric multi- 
qubit states. If the squeezing parameter < 1 (de- 
fined below), we give a quantitative relation between 
the squeezing parameter and the concurrence |25j | for 
even and odd states, where the concurrence is a mea- 
sure of the degree of two-qubit entanglement, and even 
(odd) states refer to those where only even (odd) exci- 
tations contribute. We further consider the multiqubit 
states dynamically generated from the initial state with 
all qubits being spin down via (i) the one-axis twisting 
Hamiltonian [l|, |2g, (ii) the one-axis twisting Hamil- 
tonian with an external transverse field [23, and (iii) 
the two-axis counter-twisting Hamiltonian [J . We prove 
that the states generated via the first Hamiltonian are 
spin squeezed if and only they are pairwise entangled. 
For the states generated via the second Hamiltonian and 
third Hamiltonian with even number of qubits, numeri- 
cal results for the squeezing parameter and concurrence 
suggest that the spin squeezing implies pairwise entan- 
glement and vice versa. 



II. SPIN SQUEEZING AND PAIRWISE 
ENTANGLEMENT 

A collection of N qubits is represented by the collective 
operators 



N 
i=l 



a e {x,y, z}. 



(1) 



where the Pauli operators for the i^^ qubit. The 

collective operators satisfy the usual angular momentum 
commutation relations. Following Kitagawa and Ucda's 
criterion of spin squeezing, we introduce the spin squeez- 
ing parameter 



2 (A^fl, 



4(A5«J' 



J TV ' 

where the subscript n±^ refers to an axis perpendicular 
to the mean spin {S), where the minimal value of the 
variance (A^)^ is obtained, J = N/2, and Sft^ — S ■n±. 
The inequality < 1 indicates that the system is spin 
squeezed. 

To find the relation between spin squeezing and quan- 
tum entanglement, we first give the following lemma: 
Lemma 1: For a symmetric separable state of N qubits, 
the correlation function {aiftj^ <8) f^jiij) > 0, where i and 
j can take any values from 1 to N as long as they are 
different, and am^ = (Ji ■ n± . 

Proof. We first note that the expectation values {ain^) 
and the correlation function {ain^^ crjnj_ ^ j ^-re in- 
dependent of indices due to the exchange symmetry. The 
symmetric separable state is given by 



scp 



k 



(3) 



with X^fcPfc ~ 1- The correlation function (cr,;,-j^ ® <^jrn_) 
over the separable state can be obtained from the two- 
qubit reduced density matrix 



Tr{i,2,...,Af}\{j,j}(Psop) = ^PkP^^'^ 



(4) 
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yielding 



So, the squeezing parameter becomes 



k 



(5) 



From Lemma 1, we immediately have 
Proposition 1: For an arbitrary symmetric multiquhit 
state, spin squeezing implies pairwise entanglement. 
Proof. Due to the exchange symmetry we may write the 
expectation value (S'|^) as 



Substituting the above equation into ^ leads to 



(6) 



N 



1 + {N -l){ain^(g>a■jflJ■ 



i7) 



The above equation shows that spin squeezing is equiv- 
alent to negative pairwise correlation {{cTin^^ ® ^jni) ^ 
0) |2^. This equivalence relation and the above lemma 
directly leads to the proposition. □ 

Having shown the close relation between spin squeezing 
and pairwise entanglement, we now proceed to give a 
quantitative relation between the squeezing parameter 
and the concurrence j^^. We consider an even (odd) 
pure or mixed state p. The even (odd) state refers to the 
state for which only Dicke states [2^ \n)j = | J, — J + n) 
with even (odd) n contribute. For examples, the pure 
even and odd states are given by 



I*) 



= c„|n)j, |*)o= Yl Cn\n)j, (8) 

odd n 



respectively. As we will see in the next section, these 
states can be dynamically generated via a large class of 
Hamiltonians, and can also be obtained as a superposi- 
tion of spin coherent states 20]. 

For the even and odd states, we immediately have the 
following property 

{Sp)^{S,Sp) = {Sf,S,)^Q, /3e{a;,y}. (9) 

Therefore, the mean spin is along the z direction. We 
assume that the mean spin satisfies {Sz) 7^ 0. 

With the mean spin along the z direction, we have 
n±^ — (cos^, sin0, 0), and thus the operator 5^^^ can be 
written as 



Sg — S ■ n± = cos 9Sx + sin 0Sy 



(10) 



: — min(S'fl) 
-'— min[(S'^ 



Sf. 
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{{Si + Si) - ^{si-siY-r{\s.,.Sy\ 



N 



N 



m) + \{sl)\l 



(11) 



where 5*^ = Sx i iSy are the ladder operators, and 
[^,5]+ = AB + BA is the anticommutator for opera- 
tors A and B. 

From Eq. Hll|l . we see that the squeezing parameter is 
determined by a sum of two expectation values (S*^) and 
(5^), and hence the calculations are greatly simplified. 
The larger the sum the deeper the spin squeezing. We 
also see that the squeezing parameter is invariant under 
rotation along the z direction, i.e., the squeezing param- 
eter for p is the same as that for e~^^^' pe^^^' . 

Since the inequality {SD < always holds, we 

obtain a lower bound for the squeezing parameter 



(12) 



From the above equation, we read that if = 0, 

then the squeezing parameter > 1, which implies a 
necessary condition for spin squeezing of even and odd 
states is K-S*^)] ^ 0. A direct consequence of this nec- 
essary condition is that the Dicke state \n)j exhibits no 
spin squeezing since K-S*^)! ~ 0. The associated squeez- 
ing parameter is given by 



e = i 



2n{N - n) 
N 



> 1. 



(13) 



However, Dicke states can be pairwise entangled 1291 even 
though they are not spin squeezed. 

Spin squeezing is related to pairwise correlations , and 
negative pairwise correlation is equivalent to spin squeez- 
ing [231 . Then, for our even and odd states, we have 
Proposition 2: A necessary and sufficient condition for 
spin squeezing of even and odd states is given by 



\u\-y= \{a,+ (^aj+)\ + 
where 



^ , - - > 0. 14 

4 4 ^ ' 



{ai+®(Tj+), y = ^ (1 - (cr,,^ (g) CTj^)) . (15) 



Proof. By considering the exchange symmetry, we have 
{Sl)^N{N^l)u, {Sl)^—-N{N-l)y. (16) 
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Substituting the above equation into Eq. we rewrite 
the squeezing parameter as 



e ^l~2{N^l){\u\~y) 



=1 - 2{N - 1) 



\{cr^+ (E> aj+)\ + 



(17) 



We see that spin squeezing is determined by the two cor- 
relation functions {aiz <E> (Jjz) and (cr, 
Eq. H17|) . we obtain the proposition. □ 

The two correlation functions {aiz <8) cr 
aj-^-) can be obtained from the reduced density matrix 
Pij =Tr{i 2....,Af}\{i,i}(p)- The reduced density matrix 
with the exchange symmetry is given by p9| 



CTj+). From 
) and {ai+ (g> 



x+ y y 
x+ y y 

\u X- X- 



x_ 

X* 
V- 



(18) 



in the standard basis { 1 00), |01), |10), |11)}. The following 
lemma on the reduced density matrix is useful for later 
discussions: 

Lemma 2: The matrix elements of pij can be determined 
by 



v± 
x± 

y 



N^-2N + A{S^^) ± A{Sz)iN - 1) 
" 4:N{N-1) ' 

{N-l){S+)±{[S+,Sz]+) 



2N{N -I) 



4A^(A^- 1) 



N{N -1) 



(19) 



Proof: The matrix elements can be represented by expec- 
tation values of Pauli spin operators of the two qubits. 
v± and x± are given by 



v± =i (1 ± 2(a,, 



(20) 



and u and y are given by Eq. H15|) . 

Due to the exchange symmetry, we have 



2{Sa) 

N 



{S- 



+ / 



{O'iaO- 



N{N- 1) ' 



N 

{a,+(Tjz) = 



{[S+,Szh) 



N{N - 1) 



(21) 



From Eqs. H20|) and H21|) . we may thus express the matrix 
elements of pi2 in terms of the expectation values of the 
collective operators. □ 

The concurrence quantifying the entanglement of a 
pair of qubits can be calculated from the reduced den- 
sity matrix. It is defined as 25] 



C — Ai — A2 — A3 — A4 



(22) 



where the quantities Ai are the square roots of the eigen- 
values in descending order of the matrix product 

912 = PMcTly Cr2y)P*12i<^ly <^2y)- (23) 

In if^ . denotes the complex conjugate of pi2. Note 
that we did not use the max function in the above def- 
inition of the concurrence [2^. Therefore, the negative 
concurrence implies no entanglement here. 

Both the squeezing parameter and the concurrence are 
determined by some correlation functions. So, they may 
be related to each other. The quantitative relation is 
given by 

Proposition 3: If S,^ < 1 (|w| > y) for even and odd 
states, then 



S,^^1-{N - 1)C. 



(24) 



Proof. For our state p, from Eq. and Lemma 2, it 
is found that x± — 0. Therefore, the reduced density 
matrix becomes 

v+ u* 
y y 
y y 
u w_ 



Pij 



(25) 



For this reduced density matrix (|25|) . the associated 
concurrence is given by p3| 

(2{\u\-y), if 2y< 

\2{y-^¥^), if2y> ^U^+\u\. ^^""^ 

If \u\ > y, we have 

2y<2\u\<\u\ + ^U^, (27) 

where we have used the fact 

v+v- > v± > 0. (28) 

Then, the concurrence (|26|l simplifies to 

C = 2(1^.1-2;). (29) 

By comparing Eqs. (|17|) and H29|) . we obtain the propo- 
sition. □ 

According to Proposition 3, we have 



C = 





N- 



if - 1 

J- > if < 1, 



(30) 



from which we read that (i) if the squeezing parameter 
= 1 (no squeezing) for even and odd states, then the 
concurrence is zero (no entanglement); and (ii) if < 1, 
there is squeezing, then we have a one-to-one relation 
between the spin squeezing and pairwise entanglement. 
However, for the case of > 1, the concurrence can be 
positive, and we cannot have C < as exemplified earlier 
by the Dicke states (Dicke states are simplest cases of 
even and odd states). Although the squeezing parame- 
ter > 1 implies C < is not valid in general, in the 
next section we will observe that for some even and odd 
states the squeezing parameter > 1 does imply C < 0, 
thereby establishing an equivalence between pairwise en- 
tanglement and spin squeezing. 
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III. HAMILTONIAN EVOLUTION 

Now we consider a class of states dynamically gener- 
ated from 1 0) J via the following Hamiltonian 



H = ^iSl + XSy + jiS.j:Sy + SyS-j;) + f {S z) 

with / being a function of Sz- When 

X = 7 - f{S^) - 

and 

f^ = X = f{Sz) = 0, 



(31) 



(32) 



(33) 



the Hamiltonian reduces to the one-axis twisting Hamil- 
tonian [ll l26j| and the two-axis countertwisting Hamilto- 
nian respectively. When 



x = 7 = o, fiSz) = nSz, 



(34) 



Hamiltonian H reduces to the one considered in Refs. [23, 
IM El El El m , namely, the one-axis twisting Hamil- 
tonian with a transverse field. The one-axis twist- 
ing Hamiltonian T] may be realized in variou s q uan- 
tum systems including quantum optical systems j2fil | , ion 
traps 1^3, quantum dots cavity quantum electro- 
magnetic dynamics |37j| , liq uid-state nuclear magnetic 
resonance (NMR) system |33 , and Bose- Einstein conden- 
sates 0, 1^ . Experimentally, it has been implemented 
to produce four-qubit maximally entangled states in an 
ion trap . 

The Hamiltonian exhibits a parity symmetry. 



(35) 



where M ~ Sz + J i9> the 'number operator' of the sys- 
tem, and (— 1)-'^ is the parity operator. In other words, 
the Hamiltonian is invariant under tt rotation about the 
z axis. The symmetry can be easily seen from the trans- 
formation 



{Sxj Sy, Sz)s 



- it: Sz 



( Sx , Sy , Sz ) 



(36) 



We assume that the initial density operator is chosen 
to be 



P(0)-I0)j(0|, 



(37) 



where |0)j = • - (Xjll), and state |1) denotes the 

ground state of a qubit. The density operator at time t 
is then formally written as 



Pit) 



iHt 



p(0)e 



•Ht 



(38) 



The parity symmetry of H in H35|) leads to the useful 
property given by Eq. @. For example. 



(Sx) ^Tr[Sxe 



-iHt 



=Tr[e^'^^^S'^e 
^--{Sx). 



p(0)e 



p(0)e"^^e*^*] 



(39) 



From another point view, the state p{t) is an even state 
since the Hamiltonian is quadratic in generators Sx and 
Sy and the initial state is an even state. Then, Eq. @ 
follows directly. Since state p{t) is an even state, we may 
apply the results in the last section. Next, we consider 
three representative model Hamiltonians for generating 
spin squeezing, which are special cases of Hamiltonian 
H. 



A. One-axis twisting Hamiltonian 

We first examine the well-established one-axis twisting 
model HlSi, 



Hi - fiS^, 



(40) 



for which we have 

Lemma 3: For the state dynamically generated from |0)j 
via the one-axis twisting Hamiltonian, we always have 

e < 1. 

Proof. From the results of Refs P, E^ , we have the fol- 
lowing expectation values {ft = 2/zt) 

(S'y) (N' + N- N{N - 1) cos^-2 fl) , 

{S^_) =1 (N^ + N + N{N - 1) cos^-2 fl) . (41) 
8 

Then, we obtain a useful relation for density operator 
p(t) at any time t, 

{Si ~ S'y) = {Si) - iVV4 = -NiN ~ l)y, (42) 

where we have used Eq. H19|) . From the above equation, 
we obtain 

= M^N-ir ^^^" - ^y^" + (iSx,SyU)') 



> 



{Sl Sy)"^ 



(43) 



which implies |u| > y at any time (note that ?/ > 0). 
Therefore, the squeezing parameter always satisfies ^ < 
1. □ 

Then, from Proposition 3 and Lemma 3, we obtain 
Proposition 4: For the state dynamically generated 
from |0) J via the one- axis twisting Hamiltonian, it is spin 
squeezed if and only if it is pairwise entangled. Hence 
spin squeezing and pairwise entanglement are equivalent 
for such state. 

At times for which C = 0, the state vector is either a 
product state or an A''-partite {N > 3) maximally entan- 
gled state [35I E3 which has no pairwise entanglement, 
and thus no spin squeezing. 
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FIG. 1: The spin squeezing parameter and the concurrence 
against time t for six qubits. The parameter 7 is chose to be 
1. 



B. One-axis twisting Hamiltonian with a 
transverse field 

We consider the one-axis twisting model with an ex- 
ternal transverse field described by the Hamiltonian 

H2=^iSl + nS,, (44) 

where > is the strength of the transverse field. In 
general, this model cannot be solved analytically. Nu- 
merical results show that the squeezing parameter < 1 
for the dynamically generated state exp(— ii/2t)|0) j [2^ . 
We perform numerical calculations for N from 2 to 100 
qubits, different values of fl, and /j, = 1, which indeed 
display the inequality S,^ < I- Therefore, according to 
Proposition 3, these numerical results suggest that spin 
squeezing implies pairwise entanglement and vice versa 
for such states generated from |0)j via Hamitonian H2- 
In the limit of ^ 0, the result of this subsection, of 
course, reduces to that of the previous one. 

C. Two-axis counter- twisting Hamiltonian 

Finally, we examine the two-axis counter-twisting 
model described by Hamiltonian 

H,^^iSl-S^_). (45) 

For the state generated from |0)j via Hamiltonian H3, 
the squeezing parameter can be larger than 1. Numerical 



results for N from 2 to 100 and 7=1 suggest that the 
relation (|24|) holds for even N, 

r a e^h 

C=l ^ > if < 1, (46) 

i < if > 1- 
The above equation displays an equivalence relation 
between spin squeezing and pairwise entanglement for 
states generated from |0) j via Hamiltonian H3 with even 
N. The case of TV = 6 is demonstrated in Fig. 1, which 
are plots of the spin squeezing parameter and the con- 
currence against time t. We make a conjecture that the 
spin squeezing and pairwise entanglement are equivalent 
for the states generated via the one-axis twisting Hamil- 
tonian with an external transverse field for any number 
> 2 or via the two-axis counter-twisting Hamiltonian 
for any even number of qubits. 



IV. CONCLUSIONS 

In conclusion, we have shown that spin squeezing im- 
plies pairwise entanglement for arbitrary symmetric mul- 
tiqubit states. We have identified a large class of mul- 
tiqubit states, i.e., the even and odd states, for which 
the quantitative relation of the spin squeezing parameter 
and the concurrence is given. We have proved that spin 
squeezing implies pairwise entanglement and vice versa 
for the states generated from \Q)j via the one-axis twist- 
ing Hamiltonian. For the states dynamically generated 
from |0)j via the one-axis twisting Hamiltonian with a 
transverse field for any N > 2 and the two-axis counter- 
twisting Hamiltonian with any even N , numerical results 
suggest that spin squeezing implies pairwise entangle- 
ment and vice versa. As these three model Hamiltonians 
have been realized in many physical systems, the close 
relations between the spin squeezing and pairwise entan- 
glement are meaningful and help to understand quantum 
correlations in these systems. 
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